Pulse propagation in a chain of o-rings with and without precompression 
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We implement a binary collision approximation to study pulse propagation in a chain of o-rings. 
In particular, we arrive at analytic results from which the pulse velocity is obtained by simple 
quadrature. The predicted pulse velocity is compared to the velocity obtained from the far more 
resource-intensive numerical integration of the equations of motion. We study chains without pre- 
compression, chains precompressed by a constant force at the chain ends (constant precompression), 
and chains precompressed by gravity (variable precompression). The application of the binary col- 
lision approximation to precompressed chains provides an important generalization of a successful 
theory that had up to this point only been implemented to chains without precompression, that is, 
to chains in a sonic vacuum. 

PACS numbers: 46.40.Cd,43.25.-|-y,45.70.-n,05.65.-|-b 



I. INTRODUCTION 



The dynamics of pulse propagation in granular chains 
continues to be of keen interest to physicists and engi- 
neers both because of the theoretical challenges that it 
poses, and because of the prospect of related technolog- 
ical applications. Among other applications, it has been 
argued that understanding of pulse propagation in gran- 
ular media may lead to the development of new shock 
absorbers and of instruments for detection of buried ob- 
jects Most recently, a powerful new nonlinear acous- 
tic lens consisting of ordered arrays of granular chains 
has been shown to be capable of generating high-energy 
acoustic pulses that can be used in a large variety of ap- 
plications from biomedical imaging to defense systems 
and damage detection in materials 0. While these 
applications involve propagation in higher dimensional 
granular realms, a thorough understanding of the dy- 
namics in one dimension is often pursued as a way to 
lay the groundwork for the problem, or even because 
the systems are actually composed of arrays of granu- 
lar chains. Theoretically, a granular chain appears sim- 
ple, and yet its discrete non-linear character immedi- 
ately imposes severe difficulties on the development of 
analytic methods to explain and predict the behavior of 
energy pulses in even these simplest granular environ- 
ments. As a result, the available information is obtained 
mostly from numerical and occasionally from ex- 

perimental d, results. The numerical approaches are 
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notoriously resource-intensive, and the experiments are 
constrained by the particular experimental setup. 

The numerical and experimental results have been sup- 
plemented with some analytic results that have proven 
extremely helpful and accurate. The accuracy is partic- 
ularly interesting in view of the fact that the two most 
successful theories approach the problem from two es- 
sentially orthogonal points of view, both of which have 
provided excellent results. Historically, the first was 
based on a long- wavelength approximation [TTI . IT^ that 
assumes that a chain can be treated as a continuum 
augmented with first-order corrections due to the dis- 
creteness of the system. This approach was applied to 
granular chains in which the granules just touch, so that 
there are no intergranular gaps but also no initial precom- 
pression. The long-wavelength solution to this problem 
brought to light the existence of solitary waves in granu- 
lar chains. This purely nonlinear scenario was dubbed a 
"sonic vacuum" because it supports no sound waves. The 
continuum approach has been applied in the absence [iTI - 
[Tsl as well as in the presence of dissipation [ij]. The 
success of the continuum approach is particularly note- 
worthy because it presumes a pulse width that is large 
compared to the size of the granules, and yet the ob- 
served pulse as well as the solitary wave solution that 
emerges from this approach extend over only a few gran- 
ules. The continuum methodology has been applied to 
a variety of chains with different granular configurations 
and initial conditions [Tsj . but the mathematics is quite 
cumbersome. 

A second alternative approach developed more recently 
adopts an orthogonal point of view in that it is based 
on a binary collision approximation that presumes that 
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intergranular collisions involve only two granules at a 
time. This extremely-short-wavelength approach has in 
fact been shown to yield results that are even closer to 
numerical simulation results for the "canonical" case of a 
chain of spherical granules with initially no gaps and no 
precompression [1^. Not only does this method yield 
highly accurate results for that case, but it has been 
successfully applied to a number of other configurations 
that are not easily accessible to the continuum approach, 
namely, to tapered 1^, decorated 17 1, and randomly 
decorated chains [l8| and to one-dimensional granular 
gases [l^- However, it has not yet been applied to any 
kind of precompressed chain. 

In this contribution, we study pulse propagation in a 
granular chain of toroidal rings (o-rings) placed between 
rigid cylinders that act as nonlinear springs. For such 
granules, the purely repulsive force characteristic of dry 
granular materials is modified by the topological prop- 
erties of the o-rings, leading to a hard potential propor- 
tional to the seventh power of the compression in ad- 
dition to the usual softer Hertz potential proportional 
to the compression to power 5/2. We are directed to 
this particular system because it has been studied both 
numerically and experimentally Q. Here, we make use 
of the binary collision approximation to study the pulse 
velocity in such chains of o-rings, extending the approxi- 
mation not only to a new geometry but also beyond the 
case of a sonic vacuum to precompressed chains. The 
extension of the binary collision approximation to pre- 
compressed chains provides an important generalization 
of a successful theory that had up to this point only been 
implemented for chains without precompression. 

The paper is organized as follows. In Sec.|ll]we present 
the equations of motion for the granular chain of o- 
rings and we describe the binary collision approximation. 
Next, we compare the results of this approach with those 
of the numerical integration of the equations of motion 
in Sec. IIIII Finally, in Sec. IIVI we bricfiy summarize our 
results. 



and B = 50TrDE/mdP are constants that characterize the 
elastic properties of the material Q , where d and D are 
respectively the cross-section and mean diameter of the 
o-ring, and E is the Young's modulus of the o-rings. For 
comparison with experiments, we present the values for 
these constants for the tcfion o-rings used in [sj: E = 
lAQGPa, D = 7.12mm, d = l.lQmm and m = 3.276(7. 

The granules at the ends of the chain must be con- 
sidered separately. While the leftmost granule does not 
have any granule pushing it to the right, a constant force 
Fi may be applied in this direction in the case of pre- 
compression. Therefore, its equation of motion reads 



dr^ 



F^/m^A{vi-V2f/^~B{y^-V2f. (2) 



Similarly, the rightmost granule moves according to 

Ml = _FN/m+A{yN^^-yNf'^+B{yN^^-yMf. (3) 

Without precompression, Fk = for all granules, while 
for a constant precompression, Fi = — i^w = F, and 
Ffe = for all other granules. In the case of a vertical 
chain subject to gravity with labels running downward, 
Fk ~ kmg for k ~ 1,2, . . . , N — I and F/v = —Nmg, 
g being the acceleration due to gravity. Since we are 
interested in pulse propagation, initially all the granules 
are at rest except for the leftmost or top granule, which 
has an initial velocity vq. In the case of precompression, 
either by an external constant driving force or by gravity, 
the initial positions of the granules are the equilibrium 
positions. For a chain in a sonic vacuum (that is, without 
precompression), initially the granules just touch each 
other. 

We proceed by defining the scaled variables [20l | 
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II. THE MODEL 



in terms of which the equations of motion are written as 



Consider a chain of N equal o-rings, hereafter also 
called the "granules." The double power law character 
of the elastic interaction of the o-rings leads to the fol- 
lowing equation of motion for granule k in the interior of 
the chain, 



^^A[{yk-,-ykr'^{yu^Vk+if'^] 

+ B[{yk-i - Ukf - iVk - Vk+if] + Fk/m, 



(1) 



where yk is the displacement of the fc-th granule at time 
T from its equilibrium position. Granule k is subject to 
the external force F^ , and m is the mass of the cylinders 
separating the o-rings (which is much larger than the 
mass of the o-rings). The constants A = 1.257rDE/md^^^ 



xi = fi - (xi - .Ta)^/^ - b{xi - X2f, 

Xk = fk + [{xk-i - Xkf^ - {Xk - Xk+lf^] 

+ b[{xk-i - Xkf - {xk - Xk+if], (5) 
XN = -Jn + {xn-1- xn)^^'^ + b{xN-i- xn)^, 

where a dot denotes a derivative with respect to t and 
where we have defined 
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(6) 



Consequently, the parameter b measures the relative 
strength of the two power law terms vs the initial ve- 
locity, and fk plays the role of the external force. 
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A. Binary collision approximation 

The binary collision approximation is based on the as- 
sumption that the pulse propagates through the chain by 
a sequence of binary collisions. Since this is not exactly 
the case, we must also specify the moment of passage of 
the pulse from one granule to the next. We say that the 
pulse moves from granule k to granule fc + 1 when the 
velocity of the latter surpasses the velocity of the former. 
With this approximation, instead of having to take into 
account all the equations of motion at once, we can focus 
on the interaction of just two granules at a time. Since 
the granules in the original chain are initially in equi- 
librium, so are the two granules in the binary collision 
approximation. Consequently, the equations of motion 
of these two particles during the collision may be written 
as 

£k = f -{xk~xk+if'^ -h{xi^X2)\ (7) 
ik+i = - f + {xk - Xk+if^ + b{xi - X2f , (8) 

where / is the force causing the compression of the two 
granules at the beginning of the collision. In the case of a 
sonic vacuum, this force is zero. When a constant exter- 
nal compression force is applied at the ends of the chain, 
the force is positive and / equals this external force for 
each pair of granules. In the presence of a gravitational 
force it equals kQ, the scaled gravitational constant Q 
being related to the unsealed gravitational constant g by 
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We define z — Xk^Xk+i and, for simplicity of notation, 
for now we drop the k subscripts on z and /. Subtracting 
the two equations, we have 



2/ - 2z3/2 _ 2bz^ 



(10) 



which describes the motion of a fictitious particle of unit 
mass whose displacement is z and which is moving in the 
potential 

V{z) = -2fz+-z^/^ + -bz\ (11) 
5 7 

Next we make use of conservation of energy to write 



h^ + Viz) = h',+Vizo), 



(12) 



where io = 1 is the velocity of the incoming granule at 
the beginning of the collision, and zq is obtained from 
the equilibrium condition 



f 3/2 , I, 6 

f = ZQ +bzo. 



(13) 



As discussed above, the pulse is said to reside on granule 
fc until the velocities of the two granules become equal. 
At that moment i = 0, the pulse is said to move onto the 



next granule, and the compression is maximum. Hence, 
the maximum compression Zm is obtained from the en- 
ergy conservation condition Eq. ([T2| as the solution of 

I (fil^ - ^T) + \b {^Jn - 4) -4/ (z™ - zo) = 1. (14) 

Once we know the initial and maximum compression, the 
residence time, i.e., the time spent by the pulse on a given 
granule, may be obtained as 



Z 



dz 



^l + Vizo,k)-V{z) 



, (15) 



where we have made use of Eq. and reinstated the 
subscript fc. 

In summary, in order to calculate the pulse velocity 
as it passes through granule fc, we need to find the solu- 
tions Zq and Zm of Eqs. p3p and (|14p. respectively, and 
then we must numerically integrate Eq. (jlSp (since ana- 
lytic integration appears impossible) to find the residence 
time. The pulse velocity is the inverse of the residence 
time. 



Cfc = 1/Tfc 



(16) 



III. NUMERICAL RESULTS 



In this section we compare the pulse velocity predicted 
by the binary collision approximation, Eq. (|16p , with the 
results of the numerical integration of the equations of 
motion, Eqs. ([5|). We consider three cases: chains with- 
out precompression, chains with precomprcssion caused 
by a constant force at the edges of the chain, and chains 
with precompression caused by gravity. 



A. Chains without precompression 

This is the simplest case. In the absence of precom- 
pression, zq and / vanish. The fc-independent potential 
then becomes 



V{z) = lz^/' + ^bz\ 



(17) 



and the equation for the maximum compression is sim- 
plified to 

^z^J^ + hzl = l. (18) 

These simplifications lead to the residence time 

„ . dz , , 

T= I . (19) 



, _ 8 5/2 _ 4, 7 



In the limit of small and large 6, one or the other of the 
two terms in the potential may be neglected. In these 
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FIG. 1: Pulse velocity as a function of the parameter b in the 
absence of precompression. In the main figure, the plus signs 
refer to the numerical simulation results and the open squares 
to the binary collision approximation. The dotted line is the 
plot of Eq. (pot and the dashed line is the plot of Eq. (f2T|l . 
The inset shows the relative error e. 



part of the chain) of 0.455g of mass with velocities rang- 
ing from Im/s to 800m/s. In our model the pulse is 
generated by the first granule of the chain, whose mass 
is 3.276.9 in their experiment. We adjust our initial mo- 
mentum transfered to the chain to theirs by choosing the 
parameter b appropriately, in the range [10"^", 5]. At the 
lower limit the dominant potential term is the Hertzian, 
and the pulse velocity is c = 0.82 (see Fig [IJ, which 
translates to a pulse velocity in physical units of around 
250m/s. The difference between our results and the ex- 
perimental ones in this regime is about 172m/s and is 
primarily due to the difference in the impulse generating 
method. However, it is reassuring that our results are 
within a factor of 2 of the experimental results. Exper- 
imental results were not shown for the large b limit so 
we are not able to compare with our theory. Our binary 
collision approximation predicts that the pulse velocity 
in this regime should be around llOOm/s. 



limits Zm can be calculated exactly from Eq. (fT8|) . Fur- 
thermore, the integral (|19|) can then be performed ex- 
actly leading to the low-6 pulse velocity 

and the large-5 velocity 

v^r(i/7) 

In Fig. [T] we show the pulse velocity as a function of 
the parameter b. This parameter measures the relative 
weight of the two terms in the granular interaction. For 
small values of b, the interaction is almost Hertzian while 
for large values the force is dominant. This behav- 
ior is clearly illustrated by the dotted and dashed lines, 
which are the plots of Eqs. (|20p and (HH, respectively. 
The excellent agreement of the prediction of the binary 
collision approximation is evident over the entire range 
of values of 6. In the inset we show the relative error 
£ = \icb — Cn)/cb\ between the pulse velocity predicted 
by the binary collision approximation, which we call ct,, 
and the value obtained from the numerical integration 
of the equations of motion, which we call c„. The plot 
shows that e is always smaller than 3%. Furthermore, as 
b increases the error decreases because the harder poten- 
tial becomes more important and the associated pulse is 
narrower. Consequently the binary collision approxima- 
tion turns out to be increasingly more precise, that is, the 
idea that only two granules participate in each collision 
becomes increasingly more correct jlSj . 

We end our discussion of uncompressed chains with a 
comparison of our results with the experiments of Her- 
bold and Nesterenko [Ij. In their work, the pulse was 
generated by the impact of a stainless steel sphere (not 



B. Constant precompression 

When we apply a constant force at the ends of the 
chain, thus pressing the granules together, the precom- 
pression of any pair of granules is the same, zq = S, 
as is the force / on each granule. Therefore, the pulse 
again travels with a constant (fc- independent) velocity 
along the chain. Figure [5] shows that this constant pulse 
speed increases as the precompression (or, alternatively, 
the force) increases. The inset shows the relative error, 
which increases with increasing precompression. This er- 
ror increase is a reflection of the fact that as precom- 
pression increases, more than two granules become ac- 
tively engaged in any collision event. The force caused 
by granule fc -I- 2 on granule -I- 1 for a power law po- 
tential proportional to z" is (n — 1) {Sq/Sz)^ where Sq is 
the precompression between granules k + 1 and k-\-2 and 
Sz is the additional compression between them caused by 
the traveling pulse. Hence, for large Sq, this interaction is 
enhanced and the binary collision approximation is less 
effective. The inset also shows that for larger values of 
6, the relative error e is initially smaller, as discussed in 
the previous section, but as S = So + Sz increases, the 
increase in the relative error is more prominent. Fur- 
ther, the increase in e is even faster for 6 = 20 than for 
6 = 4. Nevertheless, the relative error stays below 7% 
even for large precompressions, so our approximation is 
still useful. 

Finally, we point out that the values we have chosen 
for the precompression (up to 0.9) correspond at most 
{b = 20) to 0.9mm for the case studied in 13 (about 
13% of the diameter of the o-rings). Therefore, Fig. [2] 
probably encompasses most of the experimentally feasi- 
ble cases. Above this value, the elastic limit of the gran- 
ular interaction would not hold. 
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FIG. 2: Pulse velocity as a function of the compression. In 
the main figure, the plus signs (open squares), crosses (open 
circles) and stars (triangles) refer to the binary collision ap- 
proximation for b — 10^^, 4 and 20, respectively. The inset 
shows the relative error e for the same values of b. 



FIG. 3: Pulse velocity as a function of granule number moving 
downward along the chain. In the main figure, the plus signs 
(open squares), crosses (open circles) and stars (triangles) re- 
fer to the binary collision approximation for b — lO^'', 4 and 
20, respectively. The inset shows the relative error e for the 
same values of b. 



C. Gravitational precompression 

For a vertical chain of granules, gravity causes an ever 
increasing downward compression of the chain. That is, 
for the pair of granules k and fc + 1, the scaled force is 
f = kQ, where the scaled gravitational constant Q was 
defined in Eq. ([9|) . We again rename the initial and max- 
imum compression as ZQ^k and Zrn,k to indicate explicitly 
that these quantities now vary along the chain and hence 
depend on grain number fc. These compressions are now 
respectively the roots of the equations 

k5 = 4!k+b4,k, (22) 

"ikg {z^,k - zo,k) = ^ [z'^l - ZqS) 

+ (zo% - zl^,) - 1. (23) 

Furthermore, the full potential also becomes fc- 
dependcnt, 

V{k;z) = -2kgz+-z^/^ + -bz\ (24) 
5 7 

Therefore, the residence time Eq. ([15]), and conse- 
quently the pulse velocity, both become fc-dependent as 
well. In Fig. [3] we plot the pulse velocity as a function 
of the granule index. This dependence shows yet again 
the good agreement between the binary collision approx- 
imation and the numerical integration of the equations 
of motion. For this figure, we chose Q ~ 0.024. This 
value is too large when compared with reasonable values 
for teflon o-rings. In fact, for teflon o-rings the gravita- 
tional effects are negligible. However, it is interesting to 
understand the effects of gravity since the theory is also 
applicable to other materials. The inset shows that the 
error increases along the chain. This is not surprising, 
since the precompression also increases along the chain. 



The velocity of granules along the chain increases because 
of the increasing compression, and the binary collision 
approximation fails to fully capture this increase. Nev- 
ertheless, the effect of this failure is not very strong for 
the experimental situations that have been tested, which 
tend to involve far fewer than the 90 granules that we 
have followed in our work. Even for chains of up to 90 
granules and even with the exaggerated gravitational ef- 
fects assumed here, the relative error hardly exceeds 8% 
in the worst case scenario (small b). 

For single power-law materials, it has been shown that 
for large fc the pulse velocity scales as fc[i~i/("~i)l/2 [^^ 
[2lj |. We have verified that even for our relatively short 
chains (short for the expectation of scaling behavior), 
the pulse velocity is indeed a power law of fc. For 
b = 10'^, c - fcO ", while for 5 = 20, c - fc°-36. These 
values are not far from the values 0.17 and 0.42, respec- 
tively, which correspond to the limits of a Hertzian or a 
power 7 potential. Hence, we can see that the binary col- 
lision approximation already approaches the asymptotic 
behavior for chains of 90 granules. 



IV. CONCLUSION 

The binary collision approximation has been very suc- 
cessful in predicting pulse propagation behavior in gran- 
ular chain, providing analytic results where only numeri- 
cal ones were previously available. In this paper we have 
accomplished two goals in the further application of the 
binary collision approximation to granular chains. One 
goal has been the extension of the method to chain in 
which toroidal rings (o-rings) are placed between rigid 
cylinders that act as nonlinear springs, resulting in poten- 
tials of interaction that contain two contributions rather 
than a single Hertzian one. This system is inspired by 
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the availability of experimental results with which we can 
compare our analytic outcomes. The second, perhaps 
more important, goals is the extension of the binary colli- 
sion approximation methodology to chains with precom- 
pression, that is, beyond the sonic vacuum cases consid- 
ered in our earlier work. This extension is a challenging 
test for the binary collision approximation because pre- 
compression necessarily leads to situation in which more 
than two granules participate substantially in each colli- 
sion event. When precompression is constant, the point 
of eventual failure of the binary collision approximation 
must occur when the precompression force is sufHciently 
strong. In the case of gravitational precompression, fail- 
ure must occur when the chain is sufficiently long. How- 
ever, we find that for parameters that expansively cover 
experimental regimes the binary collision approximation 
errs by relatively little. The errors are of at most a few 
percent when resulting pulse velocities are compared with 
those obtained by numerical integration of the equations 
of motion of the full granular chain. Given the differ- 



ences in the experimental and theoretical initial setups, 
the velocities predicted by the binary collision approxi- 
mation are gratifyingly close to the experimental values 
where the latter are available. We thus conclude that 
the binary collision approximation provides a powerful 
analytic method for the study of pulse propagation in 
granular chain even in the presence of precompression. 
We continue to examine the limits of applicability of this 
powerful methodology. 
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